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IDENTIFICATION & CONTROL OF SPACECRAFT 


* THE PROBLEM 

* CONTROL 

- CLASSICAL 

- MODERN 

* IDENTIFICATION 
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CLASSICAL CONTROL 


* SCALAR FEEDBACK DESIGN 


* LEAD-LAG CONTROLLER 


K(S) 


k (s-fa) 
(s^+ 25 b+b^) 


* BANDWIDTH AROUND 1 RAD/SEC 


SHUTTLE AND PAYLOAD: FROM ALPHA TO Y-ATTITUOE 



FREQUENCY <rad/'scc> 



le m CL 


SHUTTLE AND PAYLOAD: FROH alpha to y-attitude 



MODERN CONTROL 


* START WITH FULL-STATE DESIGN 

- GOAL: MINIMIZE PAYLOAD ATTITUDE ERRORS 

- ITERATE ON CONTROL PENALTY TO ACHIEVE BW OF 5 r/s 

* DESIGN FILTER TO RECOVER LQ RESPONSE 

- USE STEIN/DOYLE ROBUSTNESS RECOVERY RESULTS 

* TEST ROBUSTNESS 
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1o0 frequency 


plant. r: Singular values of K<s>*C*Cs> q»lE8 





IDENTIFICATION 


TEST 

SIGNAL 

U 


n[uuT*^ 


MAXIMUM LIKELIHOOD ESTIMATION (MLE) 



MODEL STRUCTURE 


• State Space 

n* Modes (2x2 blocks) 


OUTPUTS 


t Frequency Domain 


(M= 


I 0 1 


jl. 


(t)+ 


[uW+w(^ 




Ci 0 ; , 
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(t) 


y(s) = (GJs) lu(s)+w(s) 1 + v(s) 


r ' 

GJ[s) = V - 2 -^-^ -2 

t Parameter Vector 

“* - 2C*^co#^, b*^, c*^; i = 1, n*) 
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MLE IDENTIFICATION SETUP 


Residual Definition 

rj^ A y(kT) - y(kT) 

Likelihood Function (Negative Log) 

A N 

L(a) = Z L. (a) 
k=0 

A ** 1 r T 1 

= ^Z^ 5 [log DET S,^ r,^J 

WHERE 

a = (af. b^, i = 1, n) = UNKNOWN PARAMETERS 

Sj^ = E^{rj^r^) = predicted RESIDUAL COVARIANCE 


HE IDENTIFICATION SOLUTION 

# Parameter Estimate (Theoretical) 

a = ArG I L(»)| = PARAMETER ESTIMATE 


• Iterative Algorithms 

Gradient: Oj+j = “j - ''L(Oj) 

Newton-Rhapson: “j+2 “ “j ' [’^L(»j)] ^vL(“j) 

VL(a) ^ (a) 


WHERE 



IDENTIFICATION ACCURACY ISSUES 


• SYSTEMATIC ERRORS: £*(&) - a* 

- MODEL-ORDER MISMATCH 

- TEST SIGNAL AMPLITUDE AND SHAPING 

- SYSTEMATIC DISTURBANCES 
-SENSOR/ACTUATOR MODEL ERRORS 

• STOCHASTIC ERRORS; a -i 

- RANDOM DISTURBANCES AND SENSOR NOISE 

- TEST SIGNAL AMPLITUDE AND SHAPING 

- IDENTIFICATION TIME INTERVAL 
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STEADY-STATE IDENTIFIABILITY ANALYSIS (YARED) 

• Expected Likelihood Function 


I («) = E. 


K<">! 


= ^ [log det S+ tr(S~^S*)] 


• Residual Covariances 

S = E„ 

S. = E 


M" 


Kalman filter predicted residual covariance matrix 
= Actual residual covariance matrix 


Note: S» and !*(«) can only be computed when the true plant parameters are 

KNOWN. 
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EXPECTED MLE IDENTIFICATION SOLUTION 


• Expected Parameter Estimate (Theoretical) 

a. - E. {aj - ArG j"«N 

• Iterative Algorithms 

Gradient; 

Newton-Rhapson: °»j+i = “*j “ 


WHERE 


(a) = 
v2[*(a) 




STEADY STATE IDENTIFICATION ACCURACY 

• SYSTEMATIC ERRORS (BIASES) 

- PARAMETER ERRORS 

“bias = “ “* = 0 MODEL MISMATCH 

- INFORMATION MEASURE (YARED) 

I(a» ; a*) - I*(a*) - I*(a*) > 0 


• STOCHASTIC ERRORS 


I [-2, ■<;.>]-' E.{[% ci>] [|!t(»>]yi-<i..] 


[v2l*( a.)] 


-1 


(N+1) 


WHEN NO MODEL MISMATCH 



STOCHASTIC ERROR WITH PROCESS NOISE 



SIMPLIFIED IDENTIFICATION ACCURACY ANALYSIS 
(FOR EACH MODE) 


e[uu^] = UI 


E[w^] = WI 
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s^+2Cti)S+o)^ 



E[vv^] = RI 




a = (lu , 2Qni, b- , . . . b , c- , , , . f c^} - {b_ or c„} 

1 m' 1' p L M 


a 2^ 

U} 

= B?! 

1 

, 2.2 
(U) ) 

SNR 

(N+1) T 

“2^ 
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(2Cto)2 

SNR 

(N+1) T 
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SNR 
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cm 
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SNR 
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A r I 

p = max { 


m j £ 1 

l,m |b|^ |c I ^ 


Ic|"[b|" U 




(MEAS. NOISE) 
(PROC. NOISE) 


(N+1) T 


(N+1) T 


HdE: THIS ANALYSIS ASSUMES THAT <-T « 1, c « 1 AND APPLIES FOR EACH MODE 




NUMERICAL RESULTS 


• PROBLEM SIZE 

12 MODES (=22 - 8R.B. - 2 SHALL) 

2 INPUTS s GIMBAL ANGLES) 

2 OUTPUTS (x> Y ATTITUDES) 

60 PARAMETERS (=12 MODES x 5 PARAMETERS/HODE) 




TEST SIGNAL. NOISE STATISTICS 
SAMPLE TIME: T = 0.1 SEC 

TEST SIGNAL! U = £|000(IN-LB)^ 

PROCESS NOISE; W = W(IN-LB)2 

MEASUREMENT NOISE: R = N x lO'^^ 


• WORST-CASE RELATIVE ERRORS AT TIME (MODE 9): 


PARAMETER 

1 SFC 

lA SEC 

(L39 HRS 


0.0265 

0.007 

0.0007 

2;oi 

3.75 

1.0 

0.1 

‘’1 

13.6 

3.6 

0.36 


2 . 7 

0.7 

0.07 

=2 

11.7 

3.1 

0.31 


SUMMARY 


• CONTROL PROBLEM 

- MODERN LQ CONTROL DESIGN WITH ROBUSTNESS RECOVERY 
PRODUCES ROBUST CONTROLLERS FOR LSS 


- ACCURATE ID ALLOWS A FIVE-FOLD INCREASE IN LOOP BW 


• IDENTIFICATION PROBLEM 

- STRUCTURAL MODES MAY BE IDENTIFIED ONE AT A TIME FOR 
SMALL DAMPING 


- LSS ID W/0 KF 

-GREATLY REDUCES PARAMETER BIASES 

-GIVES ONLY MODEST INCREASE IN STOCHASTIC ERRORS 

- RELATIVE ERRORS IN PARAMETERS AFTER ID ARE SMALLER FOR 
FREQUENCY THAN FOR DAMPING OR MODE SHAPES 

• OPEN ISSUE; HOW ACCURATE MUST ID BE FOR ROBUST CONTROL DESIGN? 
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